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Abstract: While in string theory the subject of sources in imaginary time has 
received some attention, we demonstrate the power of imaginary sources by proving 
that they constitute in several field theories a complete basis for all smooth and time 
dependent source free solutions. These proofs promote the study of imaginary sources 
to a new and crucial viewpoint for understanding time dependent backgrounds. From 
our field theory examples we further propose a completeness conjecture that every 
regular solution to a field equation of motion has a corresponding imaginary source 
configuration. We define charges for spacelike sources and show their compatibility 
with the usual charge definition for timelike sources. Many new non-singular time 
dependent field theory solutions are discussed, including Wick rotations of abelian 
instantons which demonstrate a close relationship between electric-magnetic duals 
and analytic continuation. 
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1. Imaginary Sources 
1.1 General concept 

Time-dependent backgrounds are a big issue for reconciling string theory with cos- 
mology, which has been a driving force for their study in recent years. A time- 
dependent boundary state Sen described as rolling tachyon was obtained by a 
Wick rotation, and since boundary states generically act as sources for closed strings 
including gravity, this necessarily introduces @, 0, |], [5], f| a notion of sources at 
imaginary spacetime, which we call imaginary sources in short. This concept has 
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been used in various situations implicitly, for example in Wick rotation of gravity so- 
lutions, etc. In this paper, we study how the concept of imaginary sources is general 
and fundamental. 

Inserting sources in imaginary spacetime in various field theories create many 
time-dependent classical solutions. While discovering new solutions is interesting, 
the question is how much about time dependent backgrounds do we expect to learn 
from their study. As a general remark, first of all, the fact that these solutions do 
have sources under analytic continuation is already fascinating, and more importantly 
generalizes a host of previous sources such as black holes, electrons, branes, strings to 
imaginary cases. Because of the plethora of discussions of these objects, the fact that 
their imaginary versions are consistent should give rise to much interesting study. 

A priori time- dependent solutions obtained by the analytic continuation of sources 
seem to some degree artificial. However, we claim that it is opposite: sources in imag- 
inary spacetimes have full flexibility in reproducing all the regular solutions to ho- 
mogeneous equations of motion, and will give a proof of this completeness conjecture 
for several field theories. 

Once assuming that this conjecture is true for all field theories, one important 
question is how we can recover the information of the imaginary sources from given 
time-dependent solutions. A step toward answering this question will be given as a 
new definition of charges for time-dependent configurations. 

Sources at imaginary spacetime 

Let us first review the basic properties of analytic continuation and introduce the 
concept of imaginary sources. To describe a physical system we take as our starting 
point a set of real differential equations with arbitrary number of derivatives with 
respect to {xi} for i = 1, n. The solutions of these equations are real distributions 
$ and schematically we write the situation as 

D R $ R = 0. (1.1) 

This solution can be formally promoted to the entire complex plane by changing 
all instances of real coordinates {xi} by the complex extension {zi} in both the 
differential equation and the solution 

D c $c = . (1.2) 

It is then possible to restrict the differential equation to a subset of the complex plane 
of real dimensionality n. Generally this projection produces complex differential 
equations but if this is done carefully then this restriction will generate a new real 
solution 

Dp$ P = . (1.3) 

The statement of analytic continuation can then be rephrased as follows. Take a real 
differential equation, complete it in the analytic plane and then find a real restriction 



-2- 



of this complexified system so "Analytic continuation=Completion to analytic plane 
+ projection to real subset". 



complexify /* 



\ projection 



analytic 
continuation 



The general statement of our paper is that an imaginary source can be regarded 
as follows. A solution of the differential equation with imaginary source 

D(ixi,X2, ...)&(ixx, 22, •••) = Imaginary Source S localized at (ix\ = s, x 2 = 0, ...) 

is strictly speaking a solution without source for real values of x\. For real values 
of X\ the above source 5 = 0. However the (inverse) analytic continuation of this 
solution does have a source for real values 

D(xi,x 2 , x 2 , ...) = Source S localized at (x\ = s,x 2 = 0, ...) . 

Note that we do not insist that these imaginary source solutions come from Wick 
rotations of the same Euclidean or Lorentzian metrics. Following our procedure, 
which is also discussed in Ref. 0, we shall also find a way to implicitly define an 
analytic continuation of the delta function which we otherwise would have difficulty 
in handling. 

This procedure of using imaginary sources will be very useful in generating non- 
singular solutions to technically source-less differential equation. In particular these 
solutions will generically be time dependent. This comes from the simple observation 
that if a source is placed at a complex value of time, then there is a notion of 
approaching the source and receding from the source during the evolution of real 
time. 

Completeness of the imaginary sources 

In this paper we also discover the important fact that in several classes of theo- 
ries, imaginary sources serve as a complete basis for all regular and time dependent 
solutions! Analysis of imaginary sources therefore forms a new framework for un- 
derstanding time dependence. For linear systems, one may heuristically argue for 
the validity of this statement as follows. We know that solving the homogeneous 
equations of motion with inhomogeniety given by arbitrary source, S(x), can be re- 
duced to finding the Green's function solution, G(x, xq). The general solution is then 
written as an integral of the source density times the Green's function. However if 
the source is zero over the real coordinate values then apparently this method fails. 
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Yet this is not the case. Although at first sight it appears that Green's function 
techniques are not applicable in solving the homogeneous equations of motion, on 
the contrary we will argue that in fact Green's functions can be used to solve for all 
solutions. 

The point is to properly understand the statement that sources which are zero 
valued over the real coordinate values are not necessarily zero over the complex 
values. When we write down a source S(x) it can be written as an integral over delta 
functions 



However we can examine sources which are zero over the reals but with support over 
the complex extension of the coordinates. In this case we have the identity 



when C(x,Xq) is a function with no real roots. This method clearly creates all pos- 
sible real sources and so one may reason that with the proper definition of analytic 
continuation, this method exhausts all possible ways of writing zero for real coordi- 
nate values. 

It is the goal of this paper to further describe these imaginary source techniques 
and to show the statement above can be made rigorous at least in certain cases. 
Further, this raises the possibility of a completeness conjecture for other theories 
including non-Abelian gauge theories, gravity and string theory. We will conjecture, 
based on our evidence, that for every smooth solution to the inhomogeneous equations 
of motion there is a distribution of imaginary sources. This statement, when true, 
promotes the study of imaginary sources to a crucial ingredient of all general space 
and time dependent backgrounds. 

Charges of time- dependent configurations 

Once various time-dependent solutions are obtained, a natural question is how 
we can characterize those configurations. In static cases, one can specify the solutions 
by local charges. In our strategy, solutions are obtained by Wick rotation, thus it is 
plausible that some quantity like the charge can be defined also by the Wick rotation. 
This should be a new notion because it is not conserved in the usual sense but defined 
in the full spacetime. Since this new definition of "S-charge" is obtained through 
the same type of Wick rotation used for obtaining the time-dependent solutions, it 
certainly captures the sources in the imaginary spacetime, as we will see. This new 
notion may be extended to any theory including gravity, thus possibly providing a 
new concept of "ADM S-mass" for time-dependent gravity backgrounds. There are 
other definition of the S-charges in the literature J7|, f|, H, and we will compare them 
with our definition. 
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The organization of this paper is as follows. In the next subsection we start with 
a simple solution with sources in imaginary time as a demonstration and to serve as 
a basis for the following sections. In Sec. §, we state our completeness conjecture on 
imaginary sources and give a proof for Klein-Gordon theory and Maxwell electrody- 
namics in any dimensions. In Sec. [| a new definition of "charge" capturing prop- 
erties of time-dependent backgrounds is given, and its possibilities are explored. In 
Sec. f|, we present a number of time-dependent solutions generated by the imaginary 
sources, in the Klein-Gordon systems and the Maxwell electrodynamics, including 
Wick rotation of electrons, monopoles and instantons. 



1.2 A free field theory example 

To illustrate the general concept of the imaginary source, in this subsection we will 
study imaginary sources in a simple but important example. We consider the mass- 
less real Klein-Gordon system in n dimensions. Although this example seems too 
simple, we remark that this often arises as the first order approximation to a vari- 
ety of systems from gravity to condensed matter to string theory. Interesting other 
examples will be given in Sec. £|, including electromagnetic wave solutions of various 
shapes. 

The Laplace equation of motion is 



y (— 



4>(xi,x 2 , ■ ■ ■ ,x n ) = 0. 



1.6) 



We complexify the coordinate x n into z n = x n + ix , where x is introduced as an 
imaginary part of z n . If <p{ solution to the above Euclidean equation (except 

at possible singularities corresponding to sources), then simply replacing x n by z n in 
the solution as 



(j)(x u ■ ■ • Z n ) 

gives a solution of the holomorphic equation 



;i-7) 



_d_ 

dz n 



2 n-l 

+ £ 

1=1 



d_ 

dx, 



(f)(xt, ■ ■ ■ ,x n _i,2 n ) = 0. 



Then, extracting just the imaginary part of z n we get the 'Wick rotated" solution 
of the Lorentzian wave equation of motion, 



_d_ 

dx ( 



2 n-l 



d_ 

dX; 



^1, 



x n -i,ix ) = 0. 



1.9) 



In this subsection, we see explicit examples of this procedure in 1+3 dimensions and 
find interesting time-dependent solutions which are spherical waves. 
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The first example is a Wick rotation of the simplest and the most familiar so- 
lution to the Klein-Gordon equations. To obtain a time- dependent solution in 1+3 
dimensions, we start in 4 dimensional Euclidean space. In this case we can solve the 
four dimensional Laplace equation (|1.6|) by 



N c 

^ ~ 5 x i + x 2 + 4 + (x 4 - s a ) 2 ' ^' 10 ^ 

This is the harmonic function with iV sources located at x 1 = (0,0,0, s a ). The 
charge of each source is given by the constant C a . Starting from this solution, we 
immediately obtain a time-dependent source free solution in 1+3 dimensions to the 
Klein-Gordon system by applying our above procedures. The complex extension of 
the configuration ( |1.10|) is 



N 

<P = Y~ 2 5 uP (1-11) 

^ x\ + X 2 + x\ + (z 4 - S a ) 2 
And thus we obtain a complex valued solution for the Klein-Gordon (KG) system as 

N c 

^Vt: r^^2 2 2" (1-12) 

^ ( W - S a ) 2 +X\+Xl+X\ 

Since this analytic continuation in general does not give a real valued solution, 
our procedure does not seem to make sense classically. As we will show below, 
however, the complex conjugate of the above Wick rotated solution is also a solution. 
So the real part of the above solution also solves the Klein-Gordon equation thus 
giving the real valued solution we are looking for. Physically, we are stating that 
although Wick rotating one source gives an unphysical complex valued solution, the 
Wick rotating of pairs of sources into imaginary coordinates can give a real solution. 
Adding such solutions together is possible due to the linearity of the Klein-Gordon 
equations. Even in a nonlinear theory such as string theory, pairs of imaginary D- 
branes were for example discussed in Ref. @, || £| and odd numbers of sources were 
discussed in Ref. ||. 

If we look at the above general complexification procedure, then we can un- 
derstand that we may complexify x 4 by replacing it by z 4 = x 4 — ix , not by z 4 . 
Even with this replacement, the resultant KG operator is the same. Thus, from this 
procedure we obtain a new solution to the KG equation 

N C 
^ ~ (-ix - s a ) 2 + x\ + x 2 2 + x 2 3 ^' 13 ^ 

If the parameters C a , s a are real, then this is just the complex conjugate of the 
previous solution. Thus, = 0* and the real solution is + 0*. It is also possible 
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to take the imaginary part and multiply it by i to get an alternate real solution, 
— i((f) — <p*). In other words, it is also possible to take the C a 's to be imaginary. 

To be more concrete, let us consider the simplest case in which C\ = C2 = C/2 
and Si = —S2 = s which are real; this is the same as using our analytic continuation to 
find a complex solution and then adding to it its complex conjugate. (It seems that, 
after the rotation procedure, we put sources at imaginary locations (xo, xi, £2, ^3) = 
(±is, 0, 0, 0). This is in fact true although it is not naively so. We come to this point 
later.) The solution is 



C/2 C/2 



(ixo — s) 2 + x\ + x\ + x\ ' [ixo + s)' z + xf + x\ + x\ 
C{-x\ + s 2 + x\ + x\ + x\) 



[—Xq + s 2 + x\ + x\ + x 2 ) 2 + 4s 2 xl ' 



;i.i4) 



This solution has both positive and negative valued regions. For real values of C the 
second term is just the complex conjugate of the first term, thus one may understand 
that this procedure is simply adding the complex conjugate to make the entire solu- 
tion real valued. Again, more physically stated, this is the same as Wick rotating a 
pair of similarly charged sources. 

Let us consider the singularity structure of this solution. To begin with we look 
at the locus where <ft diverges. This only occurs when the denominator is zero and 
more precisely only when the two squared terms are each separately zero 

-xl + s 2 + x\ + x\ + x\ = 4s 2 xl = 0. (1.15) 

It is easy to see that when s ^ this has no real solution. Therefore we conclude 
that the solution is regular over all real coordinate values.* 

We have plotted the function in terms of xq and x\ while putting x 2 = £3 = 0, 
in Fig. [I]. It is in fact regular, and one can see how traveling waves come in, collide 
and then go out with most of the activity centered near the light cone. This is a 



*Note however that this expectation of regular sources is not always correct. As an instructive 
example, let us consider a source at a imaginary location in the Euclidean 4 dimensional Klein- 
Gordon system, 



C 



(1.16) 



(xi) 2 + (X2) 2 + (1E3) 2 + {x4 — is) 2 
Then one finds that the singularity exists at 

2:4 = 0, and (xi) 2 + {x 2 f + (x 3 ) 2 = s 2 . (1.17) 

So this supposedly alternative method of placing sources at imaginary points has real singularities 
and does not give rise to regular solutions. In fact this solution is very unusual in that asymptotically 
it can be made to carry the same charge as a usual static source. One might wonder if this kind of 
static solution has any relationship to spinning black holes with ring shaped singularities. 
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Figure 1: A 1+1 dimensional slice of the 

solution for two sources with the same Figure 2: For two opposite charges, wave 
charge is smooth and shows a 2-sphere flips over in time which runs left to right, 
shrinking and expanding. 



non-singular collision of two wave packets traveling in opposite directions. But we 
are in 1+3 dimensions, so more precisely speaking, this is a formation of a lump by 
time-dependent shrinking 2-sphere. It shrinks to form a lump at Xq = and then 
decays to a spreading out 2-sphere. 

Let us study the sources in this Wick rotation procedure. When we complexify 
the solution as ( P-.ll ), the locations of the sources are in fact also extended in the 



complex plane. For our present example with C\ = C2 = C/2 and s\ = —S2 = s, the 
location of the singularities are simply obtained by examining where the denominator 
of the complex solution vanishes, 

(ix ± s) 2 + x\ + x\ + x\ = 0. (1.18) 

The location of the singularities in this KG system are at 



x = yjxj + xl + xj±is (1-19) 

Here we make a distinction between singularities and sources. The sources are 
point-like and at xq = ±is, they are distributed in imaginary space in terms of 
the (xo, xi, X2, X3) coordinates. These sources generate the rest of the light-cone type 
singularities. Since the singularities are all at complex coordinates values, we can 
thus describe this as saying there are no real sources in the Klein-Gordon system. 
This is the reason why we obtained a source-free time-dependent solution of the 
Klein-Gordon system. 

Having only dealt with the real part of the complexified solution, clearly there 
is also another solution which we can generate from the analytic continuation pro- 
cedure, which is to take the case of two opposite charges with G\ = — C 2 — C/2i for 
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1.20) 



real C 

_ C/2i C/2i 

x\ + x\ + x\ + {ix — s) 2 x\ + x\ + x\ + (zxo + s) 2 
Csxq 

(— Xq + s 2 + + a; 2 + x\) 2 + 4s 2 Xg ' 

This solution is non-singular with the same "imaginary" singularity type structure 
as the previous case but just with different charges. A plot of the solution is shown 
in Fig. 0. 

Note that the solution does not form a lump at the origin, is zero for xq = and 
has a maximum and a minimum. This solution can be understood as the scalar field 
theory analogue of the gravity solutions generated from diholes of Ref . |5], (| . Appar- 
ently the maximum and minimum values of the scalar field here correspond to the 
creation of large curvatures and therefore the horizons in the gravity solutions. The 
field theory solutions are different however in that they do not have an interpretation 
as the creation of a lower dimensional object as in the gravity case. 

Of course it is possible that other solutions will generate divergences on the real 
Xo axis even after analytic continuation. In such a case the Wick rotated system will 
also have a real source. For example, we can take the displacement of the sources 
in imaginary time to zero s — > 0, in which case these non-singular time dependent 
solutions to the wave equation do have a singular limit. These are the scalar field 
counterparts to the singular S-branes found in Ref. f?j] where there is a spacelike 
source existing at one moment in time. We will further discuss the usefulness of 
imaginary sources in the following sections and define a new way to calculate their 
charges. 



2. Imaginary Source Completeness 
2.1 Our conjecture 

We have seen that interesting time-dependent source-free solutions to the massless 
KG equation of motion can be obtained by applying the imaginary source technique. 
Then the following question naturally arises: what kind of solutions can we obtain 
from the technique? How general are the imaginary source solutions? Our answer 
to this question is the following conjecture: 

Any regular solution to field equations in any theory has a corresponding distri- 
bution of sources in imaginary spacetime. 

This statement is quite nontrivial. In this section, we provide a proof of this con- 
jecture for wide variety of linear systems — massless/massive Klein-Gordon theories, 
and the Maxwell theory, in any dimensions. It will be shown that for every smooth 
solution of these theories, there is a corresponding imaginary source configuration. 
Surprisingly we will find that pairs of oppositely charged sources generate all the 
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time dependent solutions in these cases, including the like charged source solutions. 
Hence dipole configurations play the crucial role as the basic ingredients of imaginary 
sources as opposed to single charges. With these insights we see that the status of 
imaginary sources is therefore promoted from curiosity to a powerful new viewpoint. 

2.2 Proof for 1+0 dimensional Klein-Gordon system 

A situation where the proof of the completeness conjecture can be explicitly and 
easily provided is the lowest dimensional example, that is, 1 + dimension. First let 
us consider the massless KG system. 

The Euclidean solution with sources at x = ±s is easily written as 

— iAx (— s < x < s) 

iAs (x<-s) (2.1) 
— iAs (s < x) 

The magnitude of the source is given just by ±zA For simplicity we take the limit 
s — > oo to get just = —iAx. The Wick rotation x — > it gives 

= At (2.2) 

which is in fact a generic solution to the equation of motion in 1+0 dimensional KG 
system once a shift in real t is included. 

The same argument holds for massive case. We put sources at x — +5 in 
Euclidean space with opposite charges, and then take the s — > oo limit to get''' 

= Ae mx + Be~ mx . (2.3) 

Then the Wick rotation gives 

= Ae imt + Be~ imt (2.4) 

which is in fact a generic solution to the massive KG equation of motion. 

After giving a proof for 1 + 1 dimensions in the next subsection, in Sec. B.4| we 
will prove the completeness of the massless/massive KG theories in any dimensions. 
Although the above proof in this subsection is straightforward, it will in fact also 
prove very useful. Before proceeding for that, here we present an issue of Wick 
rotation. Note that in the above we took s — > oo limit, but let us consider what 
happens if we didn't take this limit. It seems that this limit is unnecessary because 
the change of s does not modify the solution around the origin, = —iAx. After 
the Wick rotation we look only at the imaginary axis of x. However, this sloppy 
argument contains a subtle issue which will be explained in Sec. [4.1.1| . We have to 
take the s — > oo limit for complet eness. 

^The magnitude of the corresponding sources depends on what kind of asymptotic condition we 
chose. We may require that asymptotically the solution is zero (approaching zero exponentially), 
but this doesn't fix the magnitude of the charges completely. There remains arbitrariness of the 
choice of the charge amount, but anyway we are taking s — > oo limit and so this is not relevant to 
the following argument. 
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2.3 Proof for 1+1 dimensional massless Klein-Gordon system 

For the massless KG system in 1 + 1 dimensions, the most general solutions to the 
equation of motion are known, 

<j> = f(x-t) + g(x + t) (2.5) 

where / and g are arbitrary functions. We show that for any functions, / and g, we 
have a corresponding distribution of sources in imaginary spacetime. 

Our strategy in Sec. |1.2| for dealing with imaginary sources was to start with 
a Euclidean version of the field theory and to consider harmonic functions in that 
Euclideanized space. In two Euclidean dimensions spanned by x and y, the harmonic 
function is = log(a; 2 + y 2 ) and so an array of such charges separated in the y- 
direction produces the field 

<p = J2t Hx2 + {y ~ yn)2) - (2 - 6) 

n 

Here the source term should be given consistently by 

(dl + d 2 ) = J2^(x)S(y - y„). (2.7) 

n 

Such a potential can also be interpreted as being generated by an infinitely long 
linear static source in 1+3 dimensions. Taking the analytic continuation y — > —it of 
such a configuration we get 

^Ef^H*-^) 2 ^ 2 ) , (2-8) 

n 

for which the source is given by 

J2<lnS(y)S(-it-y n ). (2.9) 

n 

Evidently this means that the source is located at imaginary time. 

For illustration, begin with the case of two sources separated by a distance s 

= [ln(x 2 + (y - s/2f) - ln(x 2 + (y + s/2f)] . (2.10) 

The normalization of the charge has been chosen for the later purpose. The Wick 
rotation gives a Lorentzian solution which is completely regular, 

= JL \\ n (-(t - is/2) 2 + x 2 ) - ln(-(t + is/2) 2 + x 2 )] (2.11) 
2ni 



which can be simplified to 



q ts 

— arctan— —— . (2.12) 

7T x 2 -t 2 + s 2 /A y ' 
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Figure 3: Imaginary part of log wave 
with two oppositely charged imaginary 
sources at t = ±i is asymptotically flat at 
infinity It is an upside down wave coming 
in from past infinity and then it flips to 
an upright wave; time runs left to right in 
the figure. 



5 

!0-10 

Figure 4: Real part of log wave with two 
same charged imaginary sources at t = =ti; 
time runs left to right in the figure. 



This solution, as shown in Fig. |3|, is related to two opposite charged sources and is 
asymptotically flat. Now because the arctangent function is multi- valued, we have 
to specify the period of the values which the function takes. For the solution to be 
smooth also in derivatives we choose the branches of the arctangent to take 

<*/, <! (t>0) 

-1 < <f>/q < (t < 0) . V ; 

Choosing branches in this fashion amounts to three different branches of arctangent 
which smoothly fit together. The solution is created to be spatially symmetric under 
x — ► —x. However it is clearly also antisymmetric in time t —>■ —t and therefore 
imaginary time so the solution must vanish at t — 0. 

We could also take the two like charged solutions instead, whose Wick rotation 
would give, up to a proportional constant, 

<p = ln(v/(x 2 -t 2 + s 2 /4) 2 + tV) . (2.14) 

Although this solution will not turn out to play an important role in the following 
discussion, this solution is plotted in Fig. f| and is new in that the potential is not 
flat at infinity. 

Let us look at the asymptotic region of the solution ( |2.12| ) in more detail. We 



can take a limit x + t — > oo while keeping x — t finite and fixed, which gives 

q s/2 . 
(p — » — arctan . (2-15) 

7T X v 
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This shows that the asymptotic region is in fact a wave packet whose width is given 
by s and whose amplitude is not decaying but a constant. In other words, the 
asymptotic solution is in fact a (localized) plane wave. Because this plane wave has 
an arbitrary real parameter s which is the width, we can generically expect that this 
S-brane solution may form a basis for a function space in 1+1 dimensions, as the 
plane waves do. Below we present a precise proof of the completeness. 

As suggested in the generic form of the solution ( |2.5| ), one can show that the 
solution ( |2.12| ) is written as a collision of two wave packets, 

q ( -s/2 s/2\ 
6 = — arctan arctan (2.16) 

71 \ X + t X-tJ ( J 

not asymptotically but exactly. Here the arctan functions are chosen to take values 
in < arctan < tt. Let us take the thin wall limit s — > 0, then the solution becomes 

<f> = q(9{x + t)-6(x-t)). (2.17) 

The corresponding sources are distributed at (t,x) = (±s, 0), or more precisely, the 
source is given by the distribution function 

lim —5{x) [8{-it - s/2) - 5{-it + s/2)] . (2.18) 

s— >0 I 



The solution in this limit ( 2~T7 ) seems to serve as a basis of the generic solution (2.5). 
but unfortunately it leaves a constraint between the amplitudes of the left and right 
moving modes. This difficulty can not be overcome by adding combinations of the 
same sources displaced relative to each other since each source suffers from the same 
left-right mode symmetry. To overcome this difficulty, we must introduce a new 
solution whose imaginary source is located at (t, x) = (0, ±is), that is, the solution 
obtained by exchange of the roles of t and x. What we get as a result of the sum of 
those four sources is the solution (see Fig. ||) 

(f) = q(26(x + t)-l). (2.19) 

This is a purely left-going wave of a step function. Interestingly, this procedure to 
get a purely left-mover introduces the new idea of sources at imaginary space. 

At this stage it is quite easy to see that this forms a basis for generating all the 
solutions (|2.5|) , because a derivative of this step function is in fact a delta function 
5(x + t). Thus, adding various imaginary source distributed as the following source 
term 

"' /( '' '^5{x - x ) (s{-i{t - t ) - s/2) - 5{-i{t - t ) + s/2) 



lim / dx dto 



+ 2q(xo,t ) ^ t _ ^ . _ ^ _ s/2) _ _ ^ + s/2) 



(2.20) 



13 



Figure 5: Addition of the original solution with imaginary time source and 
the new solution with imaginary space source gives a purely left-moving wave. 

we may reproduce the generic solution ( |2.5| ). Here the second term is given by the 
exchange of the roles of t and x, and we let the source amplitude q and q to be 
dependent also on the real spacetime (xo, to)- In fact, if we choose 

q(x , t ) = -/'Oo - t )8(x + t ) + g'(x + t )5(x - t ), (2.21) 
q(x , t ) = f'(x - t )S(x + t ) + g'{x + t Q )8(x - t ), (2.22) 

a straightforward calculation reproduces ( |2.5| ). This is a proof of the imaginary 
source completeness in the 1+1 dimensional massless Klein-Gordon system. We see 
in this non-trivial example the power and relevance of imaginary sources as being a 
possible way to understand all features of time dependence. 

2.4 Proof in any dimensions for Klein-Gordon and Maxwell systems 

Using the results of the previous subsections, we give a simple argument to show the 
completeness of imaginary sources in any dimensions of the Klein-Gordon system 
and the Maxwell system. 

Massive Klein- Gordon 

First, let us consider the massive Klein-Gordon system. Any configuration of 
the field <ft has a representation in terms of Fourier basis, 

4>{x) = J d d+l k F(k)e ik ^. (2.23) 

For this configuration to be a solution of the equation of motion, we know that the 
plane wave is necessarily on-shell, 

k 2 Q - kj - m 2 = (i= 1,2,-.. ,d). (2.24) 

Therefore once this Fourier basis is constructed in terms of the imaginary source in 
time, a proof of the completeness is given. In fact it is very easy to obtain the cor- 
responding imaginary source distribution, if one performs a Lorentz transformation 
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on the above momentum to get 



k = ±m, ki = 0. (2.25) 
This plane wave is just the same as what we obtained in the 1+0 dimensional case of 



subsection 2.2. So if we distribute the imaginary source of the previous subsection to 



also fill all the real spatial dimensions, and perform the Lorentz transformation back, 



then we certainly get the plane wave with the momentum ( 2.24j) . Completeness of 
the imaginary sources is therefore proven in massive Klein- Gordon theory. 

Note that it is not necessary to introduce imaginary space sources in the massive 
Klein-Gordon case to remove the left/right moving symmetry as was needed in the 
case of massless Klein-Gordon. This is due to the fact that one can reduce the entire 
problem of the massive Klein-Gordon case to 0+1 dimensions which does not have 
left/right moving modes. 

Massless Klein-Gordon 

For massless KG theory, we may employ a similar argument. In this case the 
on-shell condition for the momentum of the plane wave is 

k 2 -k* = (i = l,2,...,d), (2.26) 

and a proper Lorentz transformation leads to 

k = ±h, k 2 = k 3 = --- = k d = 0. (2.27) 

Having reduced the situation to the case of 1+1 dimensions, we can apply our com- 
pleteness proof from Sec. |2.3| . There we wrote down explicitly the distribution of 
imaginary sources needed to reproduce a plane wave with momentum ( |2.27| ). Ap- 
plying the inverse Lorentz transformation gives us the desired source distribution in 
arbitrary 1 + d dimensions thereby explicitly proving the completeness conjecture. 

Maxwell system 

Interesting solutions of the Maxwell system will be given in Sec. [4.3| , but here 
we shall give a proof of the completeness for this system. The proof follows once 
we have the proof for the massless Klein-Gordon system: in a Lorentz gauge all 
the equations of motion in the Maxwell system are identical to that of the massless 
Klein-Gordon system. This is enough to know that all the Maxwell regular solutions 
can have imaginary source interpretation. One difference however is that the charge 
for Klein-Gordon is a scalar while the charge for Maxwell theory is a vector/one form. 
Lorentz transformations will therefore act differently on these charges as opposed to 
the scalar Klein-Gordon charges. 



- 15 - 



t 

particle worldline 




flux detector 



Figure 6: In Minkowski space, we can 
measure the charge of a stationary parti- 
cle by surrounding it with flux detectors 
at spatial infinity. This is labeled in the 
Penrose diagram. 




Figure 7: To calculate the charge of a 
spacelike source, we measure the flux by 
going to timelike and spatial infinity. The 
lightcone singularities generate a source at 
their intersection point. 



3. Capturing Moving Charges 

In this section we discuss a new definition of S-brane charge which we apply to 
our imaginary source solutions and find that it characterizes the imaginary source 
solutions. We show that in the correct limit it reduces to the usual charge definition 
of a static charge. This new definition of S-charge is also compared with previous 
definitions 0, |, §. 

3.1 A new definition of S-charge 

In the case of Euclidean spacetimes, we surround a point by choosing another distinct 
point and then applying the rotation group preserving the distance from the source. 
This generates a sphere surrounding the source. One can take the limit where we 
are far from the source at spatial infinity and place flux detectors there to count the 
charge. This is illustrated in a Penrose diagram in Fig. |6|. 

In the case of Lorentzian spacetimes and spacelike sources, the rotation group is 
disconnected. We wish to cover the source with hyperboloids. In fact this is a gen- 
eralization of the notion of surrounding a charge in Lorentzian signature spacetimes. 
Each connected piece of the Lorentz group will generate one of three hyperboloids. 

In 3 dimensions, we have three hyperboloids: x 2 + y 2 — r 2 = a 2 , which we call 
H 2 , is a connected hyperboloid; and x 2 + y 2 — r 2 = —a 2 gives two disconnected 
hyperboloids which we call Hi,H^ as shown in Fig. [8]. If we take the limit of these 
hyperboloids to infinity, we see that we should place flux detectors at spatial and 
timelike infinity as shown in the Penrose diagram of Fig. [7] to completely enclose the 
source. Such a notion of surrounding a spacelike charge makes sense in spacetime as 
opposed to space. 




Figure 8: To surround a spacelike source we must capture flux in the time and 
spacelike directions. 



One possible problem with these choices for hyperboloids however is that the area 
of this hyperboloid is infinite. While for a sphere a fixed radius area is f s2 r 2 sin 9d9d(p = 
47rr 2 , we have for a hyperboloid j H2 r 2 cosh 9d9d<p = oo for r 7^ 0. One point 
to note then before doing the calculation is that the induced metric on the con- 
nected hyperboloid is positive while the induced metric on the two disconnect hyper- 
boloids is negative. For the connected hyperboloid, H 2 , we have x = a cosh 9 cos 0, 
y = a cosh 9 sin <p, t = a sinh 9. For the disconnected hyperboloids Hi,H 3 we have 
x = a sinh 9 cos <p, y = a sinh 9 sin 0, t = ±a cosh 9 and we also notice that to 
cover the hyperboloids once we only need 9 G [0, 00). Calculating the area inte- 
gral we have for one disconnected hyperboloid and half of the connected hyperboloid 
— J °° a 2 sinh 9d9d(p+ J °° a 2 cosh 9d9d<p = 2Tta 2 . Note the minus sign intuitively comes 
from the fact that one are is timelike and one spacelike. If we include the contribu- 
tion from the other half we get that the total area of the three hyperboloids is 47ra 2 ! 
This is exactly then an analytic continuation of a sphere's area. 

Now let us consider a source, a delta function. For the case of Euclidean 3 
dimensional space we have that integrating the flux generated by a delta function 
over a 2-sphere is j g2 • ?r 2 sin 9d9d<p = 4n. For the Wick rotated case we have upon 
integrating the Green's function over the three hyperboloids 

/ + / + / I 4 • VV 2 (cosh9 - S inh0)d0d0 = 4tt , (3.1) 

J Hi J Hi J Hi J T 

where n 2 = —t 2 + x 2 +y 2 . This verifies the consistency of the definition of the delta 
function for spacelike sources with the new definition of S-charge. 
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3.2 Charge of the nonsingular time-dependent solutions 

Let us next try to apply this new rule of charge counting for our non-singular imag- 
inary sources. The explicit calculation is simplified by calculating at large values 
from the light cone. Remember for two static charges at x — ±s that the potential 

iS 1 1 

+ . ■ (3.2) 



yj(x — s) 2 + y 2 + z 2 yjjx + s) 2 + y 2 + z 2 

Although the configuration is non-spherical, the charge of this system is easy to 
calculate for large values of r 2 = x 2 + y 2 + z 2 where the solution simplifies to 



\/ x 2 + y 2 + z 2 



+ - (3.3) 



using a simple Taylor expansion. To calculate the charge at infinite radius, we see 
that the higher order terms do not contribute and so the charge simply is Q = 2. 

We now try to similarly calculate the charge in the case time dependent solutions 
with imaginary sources. For example the Wick rotation of the above case gives 

1 1 

+ 



y/{it - s) 2 + y 2 + z 2 y/(it + s) 2 + y 2 + z 



(( 77 2 + s 2)2 + W )l/4y ^2 + s 2)2 + 4 t 2 s 2 

where rj 1 = | — t 2 + y 2 + z 2 \. For large rj 3> s, we have an expansion 

4> = - + o(±) . (3.5) 

We see that this has a S-charge Q = 2. Despite the fact that there are neither real 
sources nor singularities, this non-singular solution has a charge! In this case the 
charge is the same regardless of how close we are to the lightcone since there are 
no real sources. On the other hand if the imaginary sources have opposite charges, 
their S-charge is zero. This is clear from Fig. ^ where we see a wave coming in and 
flipping over. The contributions from the past and future light cone regions clearly 
cancel to give zero S-charge. 

3.3 Lorentzian Gauss law 

It should be emphasized that in the case of the nonsingular solutions, our notion 
of charge seems to be conserved only in an asymptotic sense. Let us explain more 
concretely the deformations of the integration surfaces. 

In the usual definition of charges in Euclidean space, we are allowed to deform 
the integration surface arbitrarily to enclose sources. This is due to the Gauss's 
theorem. A simple statement about Gauss's theorem in 1 + 2 Lorentzian spacetime 
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is as follows. Surround the origin with a cube with sides at |(±Ax , ±Ax±, ±Ax 2 ). 
For a scalar field 0, the flux through the spatial direction x 2 is 



V0 I ^ j • £ 2 AxxAxq + V0 ( ) • (-x 2 ) A Xl Ax (3.6) 



which under Taylor expansion is 

d 2 X2 <j)(0) Ax A Xl Ax 2 . (3.7) 



The only difference for the flux calculation along the timelike direction is that the 
unit vector has negative norm x ■ x = — 1 so the flux in the x timelike direction is 

-d 2 xo (f)(0) Ax A Xl Ax 2 . (3.8) 



In total the flux through the square is 



(~d 2 X0 + d 2 Xl + d 2 jm Ax A Xl Ax 2 (3.9) 
so the generalized Gauss law in Lorentzian spacetimes (// = 0, 1, 2) is simply 

J d^<p d 3 x = Jv4>-h d 2 S . (3.10) 

Instead of just restricting ourselves to using box-like integration regions, let us 
further expand on the notion of surrounding spacelike sources. A sphere surrounds a 
point particle in Euclidean space, and a sphere in spacetime also surrounds a spacelike 
source. Arbitrary deformation of the bounding surface are allowed as long as we do 
not pass through additional charges. One can imagine then a large deformation of 
the bounding surface where the edge stretches far from the source. If the limit in 
which the contour deformation is taken to infinity carefully, then the charge bounded 
must remain constant. 

For example in the case of the singular solution, there is a concrete source produc- 
ing an associated charge. For the non-singular sources there is always an associated 
charge as calculated using our prescription. However if deformations of the hyper- 
boloids to a sphere are smooth, then we take the resulting sphere to zero radius and 
the charge accordingly vanishes due to the smoothness of our solutions. Therefore 
we should understand that our prescription of using hyperbolas for the S-charge is 
not necessarily topologically equivalent to the sphere which can shrink. The integra- 
tion surfaces consist of three hyperboloids in the case of 1+2 dimensions, and they 
are disconnected, that is why we obtained nonzero S-charge for non-singular field 
configurations. 
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3.4 Comparison with previous charge definitions 

Let us discuss relations to the previous definitions of spacelike charges. A definition 
of spacelike charge for the case of 1+1 dimensional Klein- Gordon was presented as an 
example in Ref. 0. As discussed there one choice of propagator is <p(u — t + x,v — 
t — x) = 9(u)9(v) as it satisfies 

d u d v (p(u,v) = S 2 (u,v) . (3.11) 

This solution satisfies a straightforward but probably unfamiliar Gauss law equation 
which is still of the form J v d* dcfi = f gv *d(f) but where the integral is now over the 
space and time direction. First let us calculate the S-charge of this configuration using 
the definitions of Ref. @ where we enclose the spacetime origin. If we integrate over 
a rectangle in lightcone coordinates with corners at (u,v) = (±a,±b), we calculate 
the charge of this solution to be 

Q = j*d<j> = j[5{u)6{v)du - 8{v)9{u)dv\ = Je{b)5{u)du - Je{a)5{v)dv = 2 . 

If we take this configuration to be a solution to Maxwell theory in 1+2 dimensions, 
we must also impose the current conservation law d t p + d x J x = which introduces 
the spacelike current J x = —9(x)S'(t). 

Let us compare this result with our prescription for the S-charge. In this 1+1 
dimensional case, four disconnected hyperbolas are needed to bound the source. The 
propagator $ = 9(u)9(v) can be written in terms of Milne (t, x) = ±(T cosh r, T sinh r) 
and Rindler coordinates (t, x) = ±(— T sinh r, T cosh r). In the Milne region t > \x\ 
we calculate the charge to be 

I *d$ = [ \-5{Te r )9(Te- r )e r T - 5(Te~ r )9(Te r )e~ r T]dr 

J M J 

+ [-5(Te r )9(Te- r )e r + 5{Te-' r )9{Te r )e- r }dT . (3.12) 

For simplicity take the integral to be at fixed T and from r = (— oo, oo) and change 
variables to obtain 

r°° i i 

*d$ = / -5(w)9(T/w)dw - S(w)9(T/w)dw = -9{T ■ oo) + -9(T ■ oo) = 1 



(3.13) 



' M JO 

where we have used the identity 

POO 

5(x)dx = 1/2. (3.14) 

Although the delta function is sharply defined at x — 0, we define by the meaning 
of our hyperbolic integral this symmetric condition on the integration of the delta 
function. In the second Milne region t < —\x\ we find that the answer is zero Q = 0. 
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This can be seen simply since we can obtain this second region from the first by 
simply replacing T — > — T. The final charge will then depend on 6(—Too) = giving 
zero charge. 

For the Rindler region x > \t\ we find the charge integral becomes 

J *d$ = J [5{Te~ r )6{-Te r )e- r T - 5(-Te r )6(T e - r )e r T}dr 

-[5(Te- r )6{-Te r )e- r + 5{-Te r )6{Te- r )e r }dT . (3.15) 

Again integrating at fixed T from r = (oo, — oo) and using our above definition of 
integrating delta functions 

[ *d&= [ 5(w)6(-T/w)dw - 5(w)6(T/w)dw = -6(T ■ oo) = - (3.16) 
Jr Jo 2 2 

we get a charge contribution of 1/2. Similarly in the second Rindler region the 
answer is 1/2. In total, the charge of this spacelike source is Q = 2 just like when we 
integrated before. In this case we have defined the hyperbolas to bound the charge 
which results in the same charge. 

Note that the charge using the hyperbolas would take a different value if we 
used a different way of evaluating the boundary of the integrals. Here we took only 
half of the delta function as in (|3.14|) , but if we define the integral as lim e ^ /°° dw 
then the resulting S-charge is zero. With this definition of the boundary with e, the 
hyperbolas of the integral never touch the lightcones and so with this definition it 
does not enclose the origin where the source resides. 

Before ending this subsection, we evaluate the charge of the Wick rotated log- 
arithm harmonic function in two dimensions using the definition of Ref. J7|. The 
solution to the Green's function in 2 dimensional Euclidean space is known to be 
(l/47r) ln(x 2 + y 2 ) and its Wick rotation (1/471*) ln(x 2 — t 2 ) was the basis of our dis- 
cussion in the previous section. Although we have referred to this kind of solution 
as a delta function for the wave equation due to analytic continuation, this kind of 
procedure is a bit subtle in the presence of delta functions. Let us therefore verify 
the form of singularity by applying Gauss's theorem as before. To ensure that Wick 
rotation is consistently working we check the charge of this Wick rotation using once 
again light cone coordinates (u,v). Taking into account the ie prescription which 
rewrites integrals in terms of the principal value and delta functions we have 

A W 2 .2\ Vdu - Udv in i ^\ 

*d \mx — t ) = (3.17) 



vdu I P I — ) — in5(uv) J — udv ( P ( — I — m8{uv 



Q . . . , 

UV J J \ \uv 

-in Sign(v) 5 (u) du + in Siga(u)5(v)dv] = 4ni (3.18) 

so the solution = (i/4n)ln(x 2 — t 2 ) has a unit of S-charge. 
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Figure 9: Surround two spacelike source using hyperbolas, and then cancel contributions 
whenever possible. For example contributions from hyperbolas 3 and 5 cancel. After 
making all possible cancellations we are left with a similar configuration as for one S-brane. 



3.5 Relation to usual static charge 

In this subsection we count more complicated configurations of charge. First take 
the case of two S-branes shifted in time as discussed further in Sec. [4.1.2| . Our rule 
is to surround both S-branes by a single large hyperboloid as shown in Fig. ^. These 
hyperbolas have cancelling contributions and we can deform the calculation into 
three hyperboloids again as shown in Figure |^. The calculation for two S-charges is 
similar to that of one S-charge. In the process we must integrate over the lightcone 
discontinuities. 

A further interesting observation which shows the consistency of our definitions 
is the following. Allow for a linear distribution of singular spacelike sources in the 
time direction as we will discuss in Sec. |4.1.2| . Generalizing the definition of S-charge 
just outlined above for two charges then we should integrate over three disconnected 
regions. If the S-brane is long lived then two of the regions are past the singularities 
and do not contribute much to the total S-charge. We are just left with then a 
circular integral of the usual charge for a particle, times an integral over time as 
shown in Fig. [TO 

^/S-charge Qstatic charge ' • (3.19) 

Our definition of S-charge is then related to the usual charge of a static parti- 
cle charge in the same way as static charge is related to a linear charge density in 
space, Q = XL. In this way our definition of S-charge reproduces the usual static 
charge. Our S-charge is therefore also consistent with the fact that integrating singu- 
lar spacelike sources over time reproduces the potential of an static charged particle 
as we will mention in Sec. |4.1.2| . Because spacelike sources can be thought of con- 
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cretely as the individual bits of a particle's worldline trajectory, there is something 
quite fundamental about spacelike sources in time dependent systems as opposed to 
particles. 

3.6 Secondary S-charge 



time 




In addition to the above definition of S-charge, 
there has also been a definition of S-brane charge 
proposed in Ref . [|J] . First let us review this def- 
inition. Turning away from the Klein-Gordon 
equations, consider the a positive charge and 
a negative charge forming a dipole. There are 
three charges we can calculate which are the 
dipole charge and also each individual charge. 
Place a negative/positive charge at z = ±s so 
the potential is q 



y i /x 2 +y 2 + (z—s) 2 W x 2 +y 2 +(z+s)< 

Figure 10: The integral over a small We surround the negative charge by a two sphere 

region of a hyperbola turns into an which we deform until it fills out the half plane 

integral over a cylinder with the circle between the two charges and goes out to infinity 

in the space directions and the length in the upper half plane. The electric field along 

of the cylinder along time. this z = plane is 



E, 



—c(z — s) 



+ 



c(z + s) 



2cs 



[x 2 + y 2 + (z - s) 2 ) 3 / 2 [x 2 + y 2 + (z + s) 2 f/ 2 {x 2 + y 2 + s 2 f/ 2 
If we change variables so p 2 = x 2 + y 2 then the flux through the half plane is 



.(3.20) 



E z pdpd<f) = 2tt 



2cs 



( p 2 + s 2)3/2 



pdp = 2ti{2cs) 



'p 2 + s 2 



-Anc (3.21) 



which is the same as the calculation using the usual sphere. Such a charge definition 
works only on two opposite charges and is zero when the charges are the same. This 
method picks up one of the charges but only in a dipole configuration. The flux 
calculation is of course the same for any choice of infinite half plane between the 
two charges; the only difference is that there will be also be contributions from the 
electric field in the other (x, y) spatial directions. 

Wick rotating the two oppositely charged sources by z —> ir does not affect 
the above charge calculation. What we have though is a different interpretation of 
the charge. The potential = A — > and the original electric field along the z 
direction becomes an electric field along the x A = it direction so overall F tz = F X ± T . 
Specifically the electric field along x 4 at time r = is 



[x" 



2cs 

y 2 + s 2 ) 3 / 2 



(3.22) 
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and the rest of the calculation is the same. Calculating the flux at a different time 
gives the same answer but the difference is that now there is also the presence of time 
dependent magnetic fields. If we were to calculate the charge of such a potential in a 
scalar field theory, we would obtain zero S-charge according to our prescription since 
this is a wave which flips its amplitude. 

Summarizing the difference between the definitions of S-charges, we find that our 
new definition of S-charge seems consistent with the notion of calculating the flux 
from putting in a single delta function sources in time. It is also able to calculate the 
total amount of charge which is focused into a small region of space and time. The 
notion of S-brane charge in Ref. [|J on the other hand is a way to calculate the charge 
of dipole type configurations. These two charge definitions produce different results 
in the cases we have studied. However it is very natural for spacelike sources to 
have various charges. We are well familiar with the fact that distributions of electric 
charges have monopole, dipole, quadropole, and higher multipole moments and that 
is something similar to the situation here. 

4. Explicit Examples 

In this section we provide explicit examples of the imaginary source solutions for 
massless/massive Klein-Gordon systems and also for Maxwell system. In particular, 
Wick rotations of electrons, monopoles and instantons are shown to give cylindri- 
cal/spherical electromagnetic waves. 

4.1 Arrays of massless Klein-Gordon solutions 
4.1.1 Infinite imaginary arrays 

Arrays in 4 dimensions 

In Euclidean space, let us take an infinite array of alternating point charges. The 
potential is 

; _ oo (x 4 - s/2 + 2ns) 2 + r 2 + (x 4 + s/2 + 2ns) 2 + r 2 ^ ^ 

where r 2 = (xi) 2 + (a^) 2 + i x 3) 2 - We have been dealing with finite numbered arrays 
but what about infinite numbered arrays like the one above? Does it even give a well 
defined potential when we are away from the point charges? Although an infinite 
array may seem like a very special configuration, it is important to discuss this case 
especially in relationship to the rolling tachyon. 

Let us check that 4> a rray is bounded when away from a single charge, for the case 
of real C\ = C<i- First it clear that the denominators are squares of positive numbers 
and so will be minimized when each positive number is as small as possible. This 
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means that we should take r = 0. Due to the periodicity of the solution we know that 
we can maximize the solution by simply maximizing it in the range — s/2 < £4 < s/2. 
In this range all terms with n ^ are maximized at £4 = ±s/2. Summing all such 
non-singular terms the sum is bounded by 

7 < Jr^ (2^F + (x 4 - s/2) 2 + (x 4 + s/2) 2 ' (4 ' 2) 

The infinite sum is bounded and the remaining two terms are only divergent at the 
singularity By taking absolute values, we can likewise show that arrays of positive 
and negative charges (C\ = — C2) give a bounded sum for the potential except for 
when we are at a source. 

In the case of time dependent solutions, we can similarly bound the potential. 
One easy way to check this statement is to realize that the potential for two positive 
charges was seen to maximized at the origin of spacetime (see Fig. [l]). This means 
that for each pair of same charges in the infinite sum the same is true that the 
potential is maximized at xq = 0, x\ = X2 = 23 = 0. This again reduces the sum to 
an infinite series of inverse squares which is finite. The difference is that in this case, 
all the terms are non-singular since we can never approach any of the singularities 
at xq = {imaginary values}. We conclude that the potential is completely bounded. 
As pointed out in Ref. [0], this infinite summation can be done explicitly for C\ = 
C2 = C, using the relationship 

n 1 iz 1 
tanhz = lim — -. for 4. Z (4.3) 

n->oo 27r(i — k) + Z IT 2 

k=—n v z ' 

we find the non-singular form 

C C 



•'RT 



2 



Xq — is/2 + 2ms) 2 + r 2 — (xq + is/2 + 2ins) 2 + r 

= [ tann (-(-^0 + r)j + tanh ^-(x + r) j (4.4) 

At time Xq = this function is monotonically decreasing as a function of the radius 
r from the value 1 to zero. This can be interpreted as a localized lump. For large 
values of the radius r 3> (xq,s), we have that this is a nearly static potential. 
Furthermore, this has a characteristic behavior 0kt ~ 1/r as opposed to the original 
1/r 2 contribution from each source. In the far past and far future, this potential goes 
to zero. 

Here we see that arrays of imaginary sources of codimension four, can create time 
dependent versions of real sources of codimension three. In effect we are creating a 
nonsingular time dependent version of a static point source! This can be understood 
from the fact that an infinite summation of sources is effectively an integral and 
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so the 1/r 2 behavior sums to 1/r. Imaginary source summation is not effectively 
a summation over time but a summation over space because of the "z" imaginary 
factor. Therefore essentially any infinitely extended in imaginary space summation 
of one type of charge will lead to the same kind of asymptotic 1/r behavior at spatial 
infinity. * 

A peculiar fact about such infinitely 
extended solutions is that at fixed time 
and large radius, these solutions have 
finite static charge. By "static charge" 
here we mean integrating the three form 
*d(j) over the time direction and angular 
directions at a fixed radius, f. . Q *d(f>. 
Next divide the integral by the period 
of time we integrate over and this av- 
eraged result is what we mean by the 
static charge. If we take the limit of the 
radius to infinity first, then the static 
charge is a finite constant. Therefore 
an observer would expect to see an elec- 
tron at the origin. However over a long 
enough period the time dependent terms 
have a dramatic effect and are not ig- 
norable in general. The effect of time 
dependence is to transport the field configuration to and from infinity. 

While we discussed infinite arrays of one type of charge, we can also discuss 
what happens to alternating charge arrays. In fact only these dihole charge type 
solutions have known gravity solution duals @, It is actually a simple process to 
produce the infinite array of such solutions, since the sum of alternating charges can 
equivalently be thought of as the sum of two infinite arrays of like charges. Simply 
take the above solution and translate Xq — > Xq + is/ 4 and take the imaginary part of 
the result. Using the standard hyperbolic tangent addition rule we get@ 



Figure 11: An infinite array of like codi- 
mension four charges is asymptotically in the 
spacelike direction a codimension three parti- 
cle and a decaying wave in the past and future. 
The lightcone clearly separates the character- 
istics in the two different regions. 
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cosh(27r(r + x )/s) cosh(27r(r — x )/s) 

which is clearly going to zero at large radial values and in the far past and future. 

One can also modify the charges G\ and C<i so that they are dependent on n. 
As long as the modification is mild (the charge does not grow as n), the final sum is 
convergent and well-defined. This variation corresponds to various fine structure of 
the spherical shell shrinking and e xpanding. 

*We note in that gravity, the imaginary source solution corresponding to diholcs |^], as opposed 
to an infinite number of solutions in this KG field theory case, was sufficient to generate a conical 
defect. 
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Arrays in 1 dimension 

Each source localized in one dimension generates a linear potential, = \x — a\, 
which has a derivative discontinuity. Although taking an array of such sources is a bit 
subtle due to the infinities involved, a consistent representation is the triangle wave 
function. Placing negatively charged sources at x — s/2 + 2ns and positively charged 
sources at x = — s/2 + 2ns for all integers n, then the potential is — 2x/s for x = 
[0,s/2], = 2-2a;/s for x = [s/2,3s/2] and = -4 + 2x/s for x = [3s/2,2s], and so 
on. This function periodically repeats with period 2s. The derivative discontinuities 
occur at the sources of the charge and alternate sign. This has a well known Fourier 
representation 

= — E — sin — • 46 

n 2 ^-^ n z \ s / 

n odd 

We can then easily shift x and Wick rotate to get what appears to be an infinite sum 
of Sen's rolling tachyons 

= ~2 E C ° Sh ( — ) • ( 4 - 7 ) 

n odd 



Plotting this solution as a function of t and for finite n it looks approximately linear 
near the origin before exhibiting exponential growth. The exponential growth will 
begin at earlier times as we add more and more terms so in the limit of an infinite 
number of terms there is no linear region. In fact we can easily show that this 
infinite series has a zero radius of convergence due to the hyperbolic terms growing 
too rapidly as compared to the 1/n 2 suppression factors in the coefficients. § 

One might note a particular feature of this analysis in that if we extrapolate the 
Wick rotation of the Fourier expansion for all periodic arrays of solutions, then the 
Wick rotated solution is apparently badly divergent unless there is a very special 
cancellation of the exponentially growing modes. So while 

f(x) = J2^e 2mnx (4.8) 

n 

is a reasonable expansion for periodic static functions, the Wick rotation of this func- 
tion is not necessarily well represented by the Wick rotation of the Fourier expansion 
in the case of infinite sums being needed. If we Wick rotate x — > it, then we have 
the expansion in terms of exponentially growing modes 

/(^) = E a « e " w ( 4 - 9 ) 



§ Shifting and analytically continuing x results in an infinite sum of hyperbolic sines which does 
not help the convergence properties. Also taking the sources to all have the same charge does not 
help and in fact the potential does not appear well defined in this case. 
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Figure 12: Two (left) / three (right) non-singular pairs of S-branes parallel in time. 



which is certainly not absolutely convergent and in fact leads to divergences. However 
for specially chosen coefficients a n , this sum does converge for a certain radius: an 
example is a n = e~ n2 . 

As we already know Wick rotations of infinite arrays such as the one in 4 dimen- 
sions we saw are finite. In the present case analytic continuation of the Fourier sum 
is not reliable and does not provide us a proper representation of the Wick rotated 
function. The Wick rotations of a solution and a representation of that solution 
in some basis are not necessarily consistent. In showing the completeness in 1 + 
dimension, we have been careful on this and took s —>■ oo limit to avoid the confusion. 

4.1.2 Arrays in real time 

In addition to arranging pairs of singularities at one imaginary time, we can also 
array them in time. This is very simple since for example whenever we encounter 
terms of the form — (t — is) 2 + (xi) 2 + (x 2 ) 2 + (x 3 ) 2 , we can also simply translate this 
in real time by t — > t + t . The case of two and three such configurations are shown 
below in Fig. |I2| . 

Next we ask what happens if we arrange S-branes continuously in time. If we 
arrange the singular S-branes sources in time then we can create the usual static 
objects in this way. This is simple to see from the KG example. If we have the 
equation 

(-d 2 + d\ + d\ + dl)<j) = 6(t - t )S(x 1 )5(x 2 )S(x 3 ) (4.10) 

and integrate both sides with respect to all values of to, then the right hand side 
becomes the three dimensional spatial delta function. This will also automatically 
produce the usual harmonic solution $ = Jdt <\> = l/Airy/W to the Laplacian equa- 
tion in three dimensions. In other words integrating an S-brane in time produces a 
pointlike static solution. It is also possible to create sources that begin and end in 
time, which are singular S-branes existing for finite time. Like the S-branes which 
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Figure 13: Left: Integrating the infinite array for a short period of time, —5 < t < 5, 
gives a solution similar to the un-integrated version. Right: Integrating the infinity array 
for — 10 < t < 10, we find that the solution develops a nearly constant region near the 
origin of spacetime. 



exist for a moment in time, finite sources create light cone discontinuities starting 
and ending on the endpoints of the S-branes. 

Having discussed the singular case, what static solutions we produce upon inte- 
grating imaginary sources in time? Take the four dimensional case where we must 
integrate 

dto 



(4.11) 



The result is zero, for both of our four dimensional non-singular imaginary source 
solutions ( |1.14|) and ( |1.20|) . This is technically because, for example for the latter 
solution, it takes both positive and negative values and integrating it over time gives 
zero. In addition this integral must be zero since otherwise we would have a non- 
singular static localized solitonic solution in four dimensions which is known to be 
impossible. Our result is therefore consistent. However as long as the solution is not 
infinitely extended in time (that is, the integral limits in ([4.11 ) for the time-dependent 
solution are finite) then we generate an arbitrarily long-lived solution. 

Another useful example is the partial integration of the "rolling tachyon" solution 
which we take without loss of generality from —t\ to t\ 



dt 4>(r, t - t ) 



ICtt 
r 2s V 



dto I tanh 



7r(r ~ t + t ° ) +tanh^±i 



C i 



cosh[27r(r + t x )/s\ + cosh[27rt/s] 
cosh[27r(r — t\)/s\ + cosh [27rt /s] 



(4.12) 



The spatial asymptotics are still like the un-integrated infinite array. As the 
array is integrated over longer periods of time, the potential near the origin becomes 
a constant as shown in Fig. |i~3|. 
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4.2 Massive Klein-Gordon system 

For a more explicit discussion, we here provide imaginary source solutions to massive 
Klein-Gordon theories in 1+1 and 1+2 dimensions. 

For the Euclidean massive Klein-Gordon operator, we know the Green's function 
in 3 dimensions 

-mr 

(4 - 13) 

is a solution of the equation of motion with a source, 

(d 2 -m 2 )<P = q5*(x l ), (< = 1,2,3), (4.14) 

where m is the mass and r is a radial variable in the 3 dimensional space. This 
can be used for getting a time-dependent solution in 1+2 dimensions. We can go 
ahead for getting imaginary source solutions by shifting the location of the source 
and performing the Wick rotation, 



= exp[-mA/ (-it + s) 2 + x 2 + y 2 } 
v /(-it + s) 2 + x 2 + y 2 

For example, if we take the real part of the above solution, or in other words, we 
consider like charged sources located at t — ±s, then we get a time-dependent 
solution 



[A cos(mB) - B sm(mB)} (4.16) 



r A 2 + B 2 

where A (or B) is the real (or the imaginary) part of a/ (—it + s) 2 + x 2 + y 2 , 



(A + VA 2 + 4s 2 t 2 )/2, B = -st/A, A = -t 2 + x 2 + y 2 + s 2 . (4.17) 

When an observer is spatially far from the light cone, i.e. when —t 2 + x 2 + y 2 ^> s 2 , 
then we may approximate A ~ a/— t 2 + x 2 + y 2 , so the solution looks like ~ 
gexp[— my'- 1 2 + x 2 + y 2 ]/ a/— t 2 + x 2 + y 2 . On the other hand, when x 2 + y 2 
t 2 — s 2 , we have A ~ s and so 5 ~ — £. Thus the solution is oscillatory, ~ 
gexp[— ms] sm[—mt](l/t). Finally, if we take s ^> (x, y, t) then the solution simplifies 
to the stationary wave cos(mt) (e ms /s). In sum, inside the light cone the solution is 
oscillating in time and thus looks almost like a plane wave with zero velocity, while 
outside of the light cone the solution is exponentially decaying asymptotically with 
length-scale given by proper distance from the lightcone, \/—t 2 + x 2 + y 2 . From the 
large s limit we can see straightforwardly that imaginary sources will generate all the 
sourceless 1+2 massive Klein-Gordon solutions as all such solutions arise by Lorentz 
transformations of the stationary wave. 
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Figure 14: Left: 1+1 KG solution for two imaginary source at t = ±i and with the mass 
set to one, m = 1. The ripples are in the forward and past light cones and create a lump at 
t = 0. Right: Solution for two oppositely charged imaginary source at t = ±i. The ripples 
are in the forward and past light cones and create a flat zero solution at t = 0. 



In the case of 1+1 dimensions, we start again by examining the Green's function 
in 2 Euclidean dimensions. It has the modified second Bessel function solution 

<f>{x, y) = -^i(-m^ 2 +t/ 2 ) (4.18) 

which is exponentially decreasing at large radius and goes to infinity at the origin. 
Next take two such sources at y — ±is and perform the same Wick rotation y — > it. 
Taking the real and imaginary parts of the solution to get 



> R = --^Re 

Z7T 

01 = -7T Im 

Z7T 



K^-my/x 2 - (t-is) 2 )] , (4.19) 
Kx(-my/x 2 -(t- is) 2 )] , (4.20) 



we find nonsingular time dependent solutions to the massive Klein Gordon equations 
of motion as shown in Fig. |14[ 

These solutions are asymptotically flat at large spatial directions and have damped 
oscillatory behavior in the past and future light cone. This oscillatory behavior is 
qualitatively different from the massless solutions with imaginary sources and remi- 
niscent of the minisuperspace approximation for S-branes in string theory discussed 
in Ref. ||. As with the massless KG case, these solutions are symmetric in space 
and either symmetric or antisymmetric in the time direction. 

In the case of 1+0 dimension, a generic solution has been presented in Sec. £272. 



4.3 Maxwell system 

In this subsection we consider the Maxwell equations, another important physical 
system. We restrict our attention to 1+3 or 4 dimensional space(time). The Eu- 
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clidean equations of motion are 

d m F mn = 0, (4.21) 

where m,n — 1,2,3,4. We extend the coordinate x^ to a complex variable z 4 = 
£4 + ixQ. As in the Klein-Gordon system, any solution to the above equations (with 
sources) can be extended to a solution of the complexified Maxwell equations by 
simply replacing x 4 in the solution by Z4. Restricting our attention to the X4 = 
subspace, we obtain the Lorentzian Maxwell equations 

d^F^ = (4.22) 

with //, v = 0, 1, 2, 3. Here note that we replaced the field A± by A by the rule 

A A = -iA , (4.23) 

which is necessary in accordance to the replacement d/dz 4 —>■ —id/dx^ which was 
found in the general observation in the previous KG system. Due to this replacement, 
the metric in (|4.22|) is (—1, 1, 1, 1) and no i appears in the equation. 

To apply this procedure to obtain a new Lorentzian solution from a Euclidean 
solution, we have a further requirement that the resultant Lorentzian solution is real. 
Note that for this purpose we are allowed to solve the first Euclidean equation Q4.21] ) 
with complex gauge fields. 

4.3.1 Wick rotation of electrons: cylindrical waves 

Let us present an explicit solution based on an electron-like solution in the Euclidean 
space. We adopt an ansatz A\ = A2 = A4 = for simplicity. Then the equations of 
motion of the Euclidean Maxwell system are just 

{dl + d\ + d 2 4 )A 3 = 0, d 3 A 3 = 0. (4.24) 

One example of the solutions is the electron-type, A 3 = C j \J x\ + x\ + x\, which 
is an analog of a static electron potential, where a singular source is located at 
X\ = X2 — X4 = 0. (The "real" electron can be obtained if we Wick rotate A 3 
to iAq.) Then the above Wick rotation procedure gives the following Lorentzian 
solution 

A 3 = ^=£===, A 1 = A 2 = A = 0. (4.25) 

V X l ' X 2 X 

This is a nontrivial time-dependent solution divergent at x = ±yj x\ + x\. This 
divergence, which is linearly aligned along x 3 , is propagating at the speed of light 
and lies on a light cone in the (xq,xi,X2) spacetime. 

We consider an imaginary-source extension of this solution to find nonsingular 
solutions. One might have already noticed that the situation is almost the same as 
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the massless KG imaginary source solution in 1+2 dimensions (note that it is not 1+3 
dimensions). After the Wick rotation for a pair of charges located at £4 = ix§ = ±s, 
we obtain a Lorentzian solution Aq = A± = A 3 = with A3 identified with in (|3.4| ). 
Again from the combination in 7] we expect a picture of colliding ripples as in the 
case of the massless KG system. Here however the solution is independent of x 3 and 
thus this is a colliding cylinder elongated along x 3 whose radius is shrinking (x < 0) 
and growing (xq > 0) with the speed of light. We have an energy concentration 
only on this single cylinder whose radius is changing with speed of light, so this is a 
cylindrical "pulse" , rather than an oscillating wave. Translational symmetry of this 
solution along x 3 is reflected in the location of the imaginary sources. They are at 
x = ±is, X\ = x 2 = and the sources are homogeneously distributed along x 3 . 

Using the gauge field configuration obtained, we may compute the electric and 
magnetic fields. The non- vanishing components are £3, Bi and B 2 . The config- 
uration looks quite similar to a supertube in type II superstring theory — but in 
the present case of course there is no D2-brane to support the cylinder. One would 
expect that the cylinder should shrink with the speed of light, which does coincide 
with what we have calculated. 

Another possibility is that, like the pair of the harmonic functions in (|3.4|), we 



could take opposite charges instead (or in other words, we could take an imaginary 
part of the complex solution for A 3 ). As in Sec. |4.1.2j , we may consider a continuous 
array of these solutions in time. If we continuously arrange them in time and thus 
integrate the solutions over time, we get a nonzero constant. This is different from 



the result of Sec. 4.1.2, but this is still consistent with the fact that they do solve 



the equations of motion. 

4.3.2 Wick rotation of monopoles 

For the case of "electrons," we implicitly made a double Wick rotation: first, in the 
usual electron solution we replaced the electrostatic potential by a Euclidean gauge 
field A 3 , and then we made the Wick rotation for imaginary sources by replacing A4 
by — iA . We perform a similar analytic continuation for the Dirac monopoles. We 
start with embedding a Dirac monopole configuration into 4 dimensional Euclidean 
space as 

Ai= ~f% A 2 = 9% 2 A 3 = A i = 0, (4.26) 
r(r + X4) r[r + £4) 



where r = a/ (xi) 2 + (x 2 ) 2 + (x 4 ) 2 . The embedding is into the 3-subspace spanned 
by (xi, x 2 , X4). Note that the Dirac string is at the negative axis of x 4 . Then what 
we do is first to shift the origin of the monopole from (0, 0, 0) to (0, 0, — s), and then 
to perform the Wick rotation 24 — > ixq. After these operations, the Dirac string is 
located at x\ = x 2 = and ix < —s. The important point here is that we shift the 
origin (the location of the source) in such a way that the Dirac string doesn't appear 
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in the real section of (x ,Xi,x 2 ) spacetime. Taking the singular limit where s = 
however does need more than one gauge patch as the singularity now exists in real 
spacetime. 

The gauge field A 4 becomes — iAq but this is vanishing from the first place. The 
nontrivial gauge fields A\ and A 2 become complex, so we add (or subtract) their 
complex conjugate to make them real. Since we succeed in hiding the Dirac string 
completely in imaginary spacetime, the resultant configuration is everywhere regular. 
For example, we obtain 

A 1 = -. ~ 9X2 , - + c.c. (4.27) 
r[r + s + ixq) 



where now r = a/ (x\) 2 + (x 2 ) 2 — (^o) 2 + s 2 + 2isxq. One can easily see that this is 
nonsingular, because neither r nor r + s + ix can vanish for nonzero s. 

Again, the obtained configuration is a cylindrical electro- magnetic wave (pulse), 
because the configuration does not depend on X3. 

4.3.3 Wick rotation of instantons: spherical waves 

Next, let us take an Abelian instanton as a starting point. First we review this 
singular solutions in 4 dimensional Euclidean space. The gauge field of the solution 
is written as 

/ 1 o\ 

-10 



(x) = b mn x n h(r), b 



1 
\ -1 0/ 



(4.28) 



where r = \/ x\ + x\ + x\ + x\. A straightforward computations show that 

F 12 = -2h(r) - {x\ + x 2 2 )ti(r)/r, F 34 = -2h(r) - (x 2 3 + x 2 4 )h'(r)/r, 
F13 = F 2 a = (xix 4 - x 2 x 3 )h'(r)/r, F 23 = F 41 = (x 2 x A - xix-^h' (r) /r. 

Thus the anti-self-duality equation F mn = —^€ mnpq F pq can be satisfied if we require 
an equation F 12 = — F 34 which is equivalent to —4h(r) = ^p-. This can be easily 
solved as 

h(r) = (4.29) 

This is the Abelian instanton solution with a singularity at r = 0. 
The Wick rotation procedure gives a solution 

A 1 = x 2 h(r), A 2 = —xih(r), A 3 = ix h(r), A = -ix 3 h(r), (4.30) 

with the same h(r) = C/r A but now with r = \fx\ + x\ + x\ — x\ . In this solution, 
A 3 and A are imaginary. To make a real solution, we may add a complex conjugate 
of the solution itself. When C is real, the result is, 

2C 2C 
A 1 = x 2 —, A 2 = - Xl —, A 3 = A = 0. (4.31) 
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On the other hand, if we take imaginary C = iC where C is real, then 



2C 2C 
A 1 = A 2 = 0, A 3 = -x — , A = x 3 — . (4.32) 

These two solutions are both time dependent but still singular at the spacetime 
origin r = 0. While the self-duality relation is broken when we added the complex 
conjugate, due to linearity the Maxwell equations of motion are still satisfied. 

Let us move the singularity of the original Euclidean solution a little bit, as 
before, to regularize the solution. We start with the multiple instanton solution 

N c 
A m = g W*n - trUS a ) {x , +x , + x 2l {x4 _ Sa)2)2 . (4.33) 

The Wick rotation procedure produces a new Lorentzian solution 

N r r N r r 



^ {xf + x\ + x\ + {ix -s a j z j z ' ~ ^ [x\ + xl + xl + (ix -s a ^- 
A _ Cgjixp - go) A ~iC a X3 

A3 -2^(^2:^2,^2,(-„_ Q \2\2> A 0~2^ 



f ± (xf +Xj +xl + (iXo-Sa) 2 ) 2 ' ^ (xl + xj+xl + tiXo-Sa) 2 ) 2 ' 

For simplicity we take a single location of a source (N — 1), and add the complex 
conjugate to make the solution real. Supposing that C is real, then we obtain a 
solution 

2Cx 2 (A 2 - As 2 xp -2C Xl {A 2 -As 2 x 2 ) 

1 ~ (A 2 + As 2 x 2 ) 2 ' 2 ~ (A 2 + As 2 x 2 ) 2 ' 

-2C S (A 2 - As 2 x 2 + 4a; 2 A 2 ) 8Csx Q x 3 A 
3 (A 2 + 4 S 2 z 2 ) 2 ' (A 2 + 4s 2 a; 2 ) 2 ' 1 ' 

where A = x\ + x\ + x\ — x\ + s 2 . It is easy to see that this is a regular configuration, 
because the denominator is always positive. 

Let us look at a physical property of this time-dependent solution. The energy 
concentration of this time-dependent solution can be seen clearly in its asymptotic 
expansion in which the parameter s in A is taken to be negligible, that is, (— x ) 2 + 
(%i) 2 + (x 2 ) 2 + (x 3 ) 2 ) 3 3> s 2 . Physical field strengths have the following asymptotic 
values, 

F . 4C(Qr 1 ) 2 + (x 2 f - (x 3 ) 2 + (x ) 2 ) 8Cx 2 x 3 



F 



((- Xo ) 2 + ( Xl ) 2 + ( X2 ) 2 + ( X3 ) 2 f ' id (_( a;o )2 + (;ri) 2 + (;r2) 2 +(;r3) 2 ) 3' 

—8CxiX 3 8Cx 2 Xq 



2: ' " (x ) 2 + (x 1 ) 2 + ( a ; 2 ) 2 + ( a ;3) 2 ) 3 ' 01 (-(a;o) 2 + (x 1 ) 2 + (a; 2 ) 2 + ( a ;3) 2 ) 3 ' 
Fo2= (-(x ) 2 + (x0 2 + S) 2 + (a;3) 2 )3' Fo3 = °- (435) 
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Figure 15: Left: the magnetic field F\2 of the solution (4.34) at a time slice a little after 
the collision xq = 0, in a 2 dimensional slice X3 = 0. Right: the magnetic field F13 at the 
same time and slice. This component vanishes at xq = 0. 



Note that all the entries have a Lorentz invariant length in their denominators. 
Therefore the solution ( [4.34 ) describes a 2-sphere contracting and expanding with 
the speed of light — the solution is specified by A and if one looks at an equal A 
surface one can recognize this geometrical picture. 

Let us see what happens at Xq = 0. In the solution, A vanishes, and resultantly 
the electric field vanishes; there is a magnetic field, locally concentrated near the 
origin x\ = £2 = ^3 = 0. We have plotted some of the gauge field strengths a little 
after the collision incident Xq = 0, in Fig. [15|. After a long while, the spherical wave 



(pulse) expands, as seen in Fig. [16. 




Figure 16: Left: The field strength F12 long time after the collision. One can see that 
the spherical wave (pulse) expands. Right: The same can be observed for the component 
F13. But the amplitudes of the waves are small, which is consistent with the conservation 
of the energy. 

One may think of C as a pure imaginary constant, before adding the complex 



conjugate. Then the resultant solution is different from (|4.34 ). This can be thought 
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of as "dual" of the solution ( |4.34|) , because at the collision incidence x = 0, the 
electric field is nonzero while the magnetic field vanishes. In these solutions there is 
a close and fascinating relationship between electric magnetic duality and analytic 
continuation. This is seen for example in Eqs. ([4.31 ) and (|4.32| ) which can be thought 
of as related by analytic continuation. From a short calculation one finds that these 
give rise to field strengths dual to each other. For example, F 03 of ( |4.32| ) is equal to 
F12 of ( |4.31|) if we identify C = —C. This is similar then to the usual case of electric 
magnetic duality however here we see that it applies to a new class of solutions 
coming from analytic continuations of the instanton. 

Of course, C can be an arbitrary complex constant, and thus there can be a 
mixture of these electric and magnetic collisions. Infinite arrays of such Abelian 
instantons can also be explicitly summed just like in the Klein-Gordon case using 
the same hyperbolic tangent formulas. The resulting solution is similar to a time 
dependent dyon. 



5. Further Generalizations 

In this paper we have developed the utility of complexified spacetime and its imag- 
inary sources. While it is unusual to think of fluctuations of fields to be related 
to anything else but physically accessible real spacetime objects, we conclude that 
complexification of spacetime is a useful viewpoint. Our main results include the 
proof that imaginary source configurations can be thought to be a new way to un- 
derstand solutions to homogeneous differential equations since we are able to prove 
in several classes of theories that every smooth solution has at least one correspond- 
ing imaginary source configuration. In mathematical terms our result can be viewed 
in the linear cases as unifying the analysis of singular and non-singular solutions to 
differential equations by extending Green's function techniques. Therefore we expect 
the methods of imaginary sources to be a useful and practical tool in many other 
settings both linear and nonlinear. The following applications might be interesting, 

• Non-Abelian gauge theory • Dirac-Born-Infeld theory 

• Gravity • String theory 

as all are nonlinear in equations and so the reality condition will be difficult to be 
satisfied. But it is possible that in special examples the reality can be ensured. For 
example in the case of gravity, this should be explicitly possible in all Weyl type 
ansatzes. In string theory, imaginary D-branes in pairs @, [|, |J or odd numbers || 
do give rise to sensible backgrounds. From our discussion it is clear that imaginary 
D-branes do not exhaust the physical closed string backgrounds even at the linear 
massless level. To produce independent dilaton and gauge field configurations, at 
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minimum it is necessary to use non-BPS imaginary D-branes. The notion of imagi- 
nary strings must also be developed to produce the source free background B-fields. 
A notion of imaginary strings is in fact probably already necessary to understand 
the dynamics of imaginary D-branes. 

To what extent are our current results related to such nonlinear theories. In 
this paper we provided proofs of the completeness of the imaginary sources even 
for the Maxwell system. One should note that the Klein-Gordon equation arises in 
many important cases such as the linearizations of gravity and string theory. To the 
extent that we trust that the linear approximations have corresponding nonlinear 
exact solutions, we have proved that imaginary sources are relevant in gravity and 
string theory. However this is not so straightforward. One might argue that in non- 
linear theories such as field theories with nontrivial interaction terms there will be 
no notion of "basis" because field configurations can not simply be added together 
to produce a new configuration. We hope though that S-charges associated to given 
time-dependent solutions (when the completeness conjecture is correct it always ex- 
ists) would be an interesting and intrinsic ways to characterize field configurations. 

In addition to understanding the realm of applicability of the imaginary source 
conjecture in other theories, there are many interesting directions to pursue. For ex- 
ample in gravitational theories one should develop the notion of ADM S-mass. One 
can also try to describe the interaction of imaginary sources with each other and 
real sources. Describing such processes would in fact lead to an alternate interpreta- 
tion of positron electron annihilation. Such unusual interactions between imaginary 
sources and real sources should arise since we have shown asymptotically they can 
have the same mass and charge so these quantum numbers can be conserved in inter- 
actions. A further natural question is can one generalize the notion of spontaneous 
symmetry breaking to the case of imaginary source solutions. This would include 
the development of solutions where the delta functions can be field dependent and 
take imaginary values in contrast to delta functions which are just spacetime de- 
pendent. Moreover, the arrays infinitely extended in time give rise to constant field 
configurations. If one could produce similar imaginary D-branes extended in time, 
they might be a useful clue to describing field theory duals to gravity theories in flat 
space. Finally it should be interesting to understand if notions such as T-duality are 
applicable to the case of imaginary sources. 
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